Boosting jet power in black hole spacetimes 
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The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mecha- 
nism has long been understood as an important component in models to explain energetic jets from 
compact astrophysical sources. Here we show more generally that the kinetic energy of the black 
hole, both rotational and translational, can be tapped, thereby producing even more luminous jets 
powered by the interaction of the black hole with its surrounding plasma. We study the resulting 
Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter 
case, we find that increasing the orbital angular momenta of the system and/or the spins of the 
individual black holes results in an enhanced Poynting flux. 



I. INTRODUCTION 

Enormously powerful events illuminate the universe 
that challenge our understanding of the cosmos. Indeed, 
intense electromagnetic emissions of order > 10 51 ergs [T]- 
[3] have been routinely observed in supernovae, gamma 
ray bursts (GRBs), and active galactic nuclei (AGNs), 
for example. However, despite important theoretical and 
observational advances, we still lack a thorough under- 
standing of these systems (e.g. [4]). Intense observational 
and theoretical efforts are ongoing in order to unravel 
these fascinating phenomena. While the full details re- 
main elusive, one of the natural ingredients in theoreti- 
cal models is the inclusion of a rotating black hole which 
serves to convert binding and rotational energy of the 
system to electromagnetic radiation in a highly efficient 
manner. The starting point for these theoretical mod- 
els can be traced back to ideas laid out by Penrose [5] 
and Blandford and Znajek (BZ) [6], which explain the 
extraction of energy from a rotating black hole. These 
seminal studies, along with subsequent work (see refer- 
ences in e.g. [ZH9]), have provided a basic understanding 
of highly energetic emissions from single black hole sys- 
tems interacting with their surroundings. 

Recent work has indicated that related systems can 
also tap kinetic energy and lead to powerful jets [TO] . 
This work concentrated on non-spinning black holes mov- 
ing through a plasma and highlighted that relative black 
hole motion alone (with respect to a stationary electro- 
magnetic field topology at far distances from it) can in- 
duce the production of jets. Furthermore, subsequent 
work [TT] demonstrated that even black holes with mis- 
aligned spins with respect to the asymptotic magnetic 
field direction induce strong emissions with power com- 
parable to the aligned case. These studies suggested that, 



independent of their inclination, astrophysical jets might 
be powered by the efficient extraction of both rotational 
and translational kinetic energy of black holes, and in- 
ducing even more powerful jets than the standard BZ 
mechanism would suggest. 

Galactic mergers provide a likely scenario for the pro- 
duction of the binary black hole systems considered 
here fl2j[T3]. I n sucn a merger, the supermassive black 
holes associated with each galaxy will ultimately form a 
binary in the merged galaxy. A variety of interactions 
will tighten the black hole binary. Eventually the dy- 
namics of the system will be governed by gravitational 
radiation reaction which drives the binary to merge. The 
circumbinary disk will likely be magnetized and thereby 
anchor magnetic field lines, some of which will traverse 
the central region containing the binary. Preliminary 
observational evidence for supermassive black hole bi- 
naries resulting from galactic mergers has already been 
presented [T4] , 

An ambient magnetic field threading a spinning black 
hole populates a low density plasma surrounding the 
black hole as explained by BZ [6 . Even for black holes 
with no spin, it was recently shown that the orbiting bi- 
nary interacting with the surrounding plasma can lead 
to a collimated Poynting flux [TO]. In this work, we 
consider this basic paradigm of energy extraction from 
black holes with the additional complexity of intrinsic 
black hole spin. Binaries consisting of spinning black 
holes demonstrate similar, albeit energetically enhanced, 
phenomenology. Furthermore, we investigate the depen- 
dence of the energy flux on the black hole velocity and 
highlight a resulting strong boost in the emitted power. 

In addition to exquisite and powerful electromagnetic 
detectors, soon gravitational waves will be added to the 
arsenal of phenomena employed to understand our cos- 
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mos. These studies suggest excellent prospects for the 
coincident detection of both electromagnetic and gravita- 
tional signals from binary black hole systems. Certainly, 
dual detection of electromagnetic and gravitational wave 
signals would transform our understanding of these sys- 
tems and lead to the refinement of theoretical models 

(e.g [MED- 

In the remainder of this paper, we elucidate the ba- 
sic phenomenology arising from the interaction of bi- 
nary black hole systems immersed within a plasma en- 
vironment. We focus in particular on understanding 
the Poynting flux emissions from such binaries as well 
as single, possibly spinning, black holes moving through 
plasma. We describe our equations and assumptions em- 
ployed together with some of the details of our numerical 



implementation in Sec. [TTJ Sec. [El] describes our results 
for both single and binary black holes and we provide 



concluding comments in Sec. IV 



II. IMPLEMENTATION DETAILS 

The combined gravitational and electromagnetic sys- 
tems that we consider consist of black hole spacetimes 
in which the black holes can be regarded as immersed 
in an external magnetic field. Such fields, as mentioned 
above, will be anchored to a disk. We consider this disk 
to be outside our computational domain but its influence 
is realized through the imposition of suitable boundary 
conditions on the corresponding electric and magnetic 
fields on the boundaries of our domain [10] [TT]. With 
regard to the magnet osphere around the black holes, we 
assume that the energy density of the magnetic field dom- 
inates over its tenuous density such that the inertia of the 
plasma can be neglected. The magnetosphere is therefore 
treated within the force- free approximation [6j [19] . We 
note that the contribution of the energy density of the 
plasma to the dynamics of the spacetime is negligible 
and we can ignore its back reaction on the spacetime. 

We use the BSSN formulation [2Q1 EQ of the Ein- 
stein equations and the force-free equations as described 
in [10] [TT] . We discretize the equations using finite differ- 
ence techniques on a regular Cartesian grid and use adap- 
tive mesh refinement (AMR) to ensure that sufficient 
resolution is available where required in an efficient man- 
ner. We use the HAD computational infrastructure, which 
provides distributed, Berger-Oliger style AMR [22J [23] 
with full sub-cycling in time, together with an improved 
treatment of artificial boundaries [24 ]. Refinement re- 
gions are determined using truncation error estimation 
provided by a shadow hierarchy [25], which adapts dy- 
namically to ensure the estimated error is bounded by 
a pre-specified tolerance. Typically our adopted values 
result in a grid hierarchy yielding a resolution such that 
40 grid points in each direction cover each black hole. 
We use a fourth order accurate spatial discretization and 
a third order accurate in time Runge-Kutta integration 
scheme [26]. We adopt a Courant parameter of A = 0.4 



so that Ate — 0.4Aa^ on each refinement level £. In tests 
performed here for the coupled system (and in our previ- 
ous works for the force- free Maxwell equations), the code 
demonstrates convergence while maintaining small con- 
straint residuals for orbiting black holes. Furthermore, 
we obtain for orbiting black hole evolutions agreemend 
with runs from other codes for the same initial data. 

To extract physical information, we monitor the 
Newman- Penrose electromagnetic ($2) and gravitational 
(^4) radiative scalar s [27 . These scalars are computed 
by contracting the Maxwell and the Weyl tensors respec- 
tively, with a suitably defined null tetrad (as discussed 
in [21]), 

$2 = F ab n a rh b , ^ 4 = C abcd n a rh b n c rh d , (1) 

and they allow us to account for the energy carried off by 
outgoing waves at infinity. The scalar $2 (essentially the 
radial component of the Poynting vector) provides a mea- 
sure of the electromagnetic radiation at large distances 
from an isolated system. However, as the system studied 
here has a ubiquitous magnetic field, special care must 
be taken to compute the energy flux. We account for this 
difficulty by subtracting the scalar $0 = —F a i)l a m b from 
$2. Hence, from here forward, by $2 we mean the differ- 
ence $2 —> $2 — The luminosities in electromagnetic 
and gravitational waves are given by the integrals of the 
fluxes 
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We assume that the black holes are immersed in an 
initially constant magnetic field, such as one produced 
by a distant disk surrounding the black hole system. For 
the binary systems we consider, the orbital plane of the 
evolution is assumed to be aligned with that of the cir- 
cumbinary disk. The magnetic field is anchored in the 
disk with its associated magnetic dipole aligned with the 
orbital angular momentum (chosen along the z direc- 
tion). Initially, the magnetic field is set to be perpen- 
dicular to the velocity of the black holes and the electric 
field is set to zero. Because the electromagnetic field is 
affected by the spacetime curvature, it will be dynami- 
cally distorted from its initial configuration and generate 
a transient burst of both gravitational and electromag- 
netic waves as it settles into a physically relevant and 
dynamical configuration. 

The initial magnitude of the magnetic field Bq is cho- 
sen to be consistent with astrophysically relevant val- 
ues [29] [30] . We present our results for a field strength 
which is bounded by the the Eddington magnetic field 
strength 5-6x 10 4 (M/10 8 M ) -1 / 2 G [31. As noted 
above, for these values the plasma's energy is several or- 
ders of magnitude smaller than that of the gravitational 
field. Thus, although the plasma is profoundly affected 
by the black holes, it has a negligible influence on their 
dynamics. 
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FIG. 1: Zoom-in of the electromagnetic emission |r<3>2| 2 from 
a boosted black hole (with p — 0.15x). The left frame shows 
the non-spinning black hole case (a = 0), while the right 
frame shows the spinning case (a = 0.6). The color scale is 
the same for both frames and the black hole is boosted to 
the right. The stronger emission in the vertically collimated 
region of the right frame is apparent. 



III. RESULTS 

Here we discuss our results for both boosted and 
binary black hole systems. For simplicity we introduce 
the notation L 43 = 10 43 ergs/s and compute dimension- 
ful quantities with respect to a representative system 
with total mass 1O 8 M immersed in a magnetic field 
with strength 10 4 G. We explicitly provide proportion- 
ality factors for calculations for other configurations. 
Quantities in geometric units are calculated by setting 
G = c = 1. 

Single black holes 

For single, boosted black holes, we adopt a computa- 
tional domain defined by [— 320 m&, 320 ra&] 3 and assume 
a fixed "bare" mass of the black hole (the mass it would 
have if it were static and in isolation) of = 1 in geo- 
metric units. We set the initial linear momentum of the 
black hole to p — 0.05 ni (with n an integer that we 
vary between and 4), and we set the intrinsic angular 
momentum parameter a = J/M 2 to be or 0.6 k. For 
these parameters we study the resulting electromagnetic 
collimated flux energy and its dependence on both boost 
velocity and black hole spin. 

Fig. [T] illustrates the qualitative features of the elec- 
tromagnetic emission for a boosted black hole with and 
without spin. A collimated emission is clearly induced 
along the asymptotic magnetic field direction. Signifi- 
cantly, an even stronger emission is obtained for the spin- 
ning case as expected. 

To achieve a more quantitative understanding, we com- 
pute a measure of the collimated electromagnetic lumi- 
nosity, ^collimated 7 by integrating the flux over a 15° cone 
of a spherical surface centered along the moving black 
hole with a radius r = 20 (roughly equivalent to 
40 Rh , where Rh is the Schwarzschild radius of the black 
hole). We present our results with respect to v (the mea- 
sured coordinate velocity of the black hole) and plot the 



collimated luminosity achieved once the system reaches 
a quasi-stationary state (in other words, after the initial 
transient stage) for the spinning and non-spinning cases 
in Fig. [2] 

Several key observations are evident from the figure. 
For v = 0, the electromagnetic energy luminosity does 
not vanish for the spinning black hole, while it does van- 
ish for the non-spinning black hole. This is expected, 
as the spinning black hole interacts with a surround- 
ing plasma and radiates by the Blandford-Znajek mech- 
anism [6 j. This luminosity results from the plasma's 
ability to extract rotational energy from the black hole 
and power a jet with an energy luminosity scaling as 
L BZ ~ tt 2 H B 2 [TOES] (with tt H = a/(2R H ) the rota- 
tion frequency associated with the black hole). 

For v ^ 0, both the spinning and non-spinning black 
holes have a non-trivial associated energy flux. In the 
latter case this flux arises solely from the ability of the 
system to tap translational kinetic energy from the black 
hole, while the former results from the extraction of both 
translational and rotational kinetic energies. 

The non-spinning black hole shows the expected v 2 
dependence in the electromagnetic energy luminosity 
(see Fig. [2| consistent with the membrane picture of 
a black hole as a conductor. Indeed, as was already 
indicated in the spinless case [TUl [TTJ [33j [34] , a black hole 
moving (with speed v) through an ambient magnetic 
field acquires an induced charge separation (oc v). The 
membrane paradigm [35] of a black hole explains this 
induced charge by regarding the black hole as a con- 
ductor moving through a magnetic field, and hence the 
induced charge is analogous to the classical Hall effect. 
With this observation and the induction equation, it is 
straightforward to conclude an electromagnetic energy 
flux will be produced with magnitude oc v 2 B 2 . Such 
a quadratic dependence on speed is apparent in the 
figure. These results are also consistent with the work 
of [36] which studied the Poynting flux associated with a 
moving conductor in a magnetized plasma, in particular 
the motion of artificial satellites in orbit. They found 
that this flux obeys L v « (v/v^^B 2 [36] where v a u, 
the propagation speed of the Alfven modes, is the speed 
of light in a force-free environment. Furthermore, a 
misalignment is expected between the collimated energy 
flux and the original magnetic field orientation such 
that tan(a) = v/v & \{. With the cautionary note that 
measuring this angle is a delicate enterprise within Gen- 
eral Relativity, especially in the strongly curved region 
around the black hole, such a relation is manifested 
in our results. For instance, for a = 0, p x = 0.10 the 
measured angle obeys tan(a) = 0.07 instead of the 
predicted 0.08 [36]. 

The luminosity for the spinning black hole also demon- 
strates the same quadratic dependence in velocity, with 
the additional component from the spin. Notice that 
the difference between the obtained luminosities in the 
spinning and spinless cases remains fairly constant for 
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the different values of v. Thus, to a reasonably good 
approximation, we can express the luminosity obtained 
as a sum of a spin dependent component and a speed 
dependent one, L C oiiimated - L spin + £ sp eed, the latter 
scaling as L spe ed = L 2 v 2 . A fit to the case a = 
gives L2 = 127 (M$ B4) 2 L43. Next, using this value for 
the case a = 0.6, we obtain the fit for L spm = L\ = 
0.87 (M 8 5 4 ) L 43 . We can then provide an estimate of 
the associated luminosity for general cases as follows. 

The spin component for the general case can be ex- 
pressed in terms of its dependence on the the black hole 
spin as L sp m = Li(f2#(a)/r2#(a = 0.6)) 2 . The boost 
component takes the already mentioned quadratic form 

^speed = L2V 2 . 

Therefore, we have for the estimate 
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Notice that for a = 0.6, the non-rotational contribution 
^speed to the emitted power becomes larger than the 
rotational one for speeds approximately v > 0.08c. This 
relationship, for example, would predict a luminosity for 
an a = 0.95, v = 0.5c black hole to be ~ 36L 43 (M 8 B A ) 2 . 
This phenomenology strongly suggests that the Poynting 
flux can tap both rotational and translational kinetic 
energies from the black hole and that faster and more 
rapidly spinning black holes have a stronger associated 
power output. 




FIG. 2: Collimated luminosity as a function of the boost ve- 
locity of the black hole for two different spin magnitudes. The 
dashed lines are obtained from a quadratic fit of the form 
^collimated = ^spin + L 2 v 2 where L sp i n and L 2 are two con- 
stants. The constant L 2 is the same for both fits, supporting 
the argument that the luminosity does contain two different 
components, one for spin and one for boost. 



Binary black holes 

We turn our attention now to the orbiting binary black 
hole case. We consider equal-mass configurations in 
which each bare black hole has mass = 0.483 and 
the pair is separated by a distance D w 16 m^. We adopt 
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J ADM 
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(0.00208,-0.11235) 





1.0 


0.108 


u/d 


(0.00208,-0.11235) 


±0.515 


1.1 


0.108 


u/u 


(0.00208,-0.11235) 


0.515 


1.1 


0.349 



TABLE I: Binary black hole configurations considered here. 
All begin with no momentum in the z direction. 



a computational domain of [— 638 m^, 638 ra^] 3 and con- 
sider configurations with individual spins either vanish- 
ing or in the up (+z) or down (—z) orientation with spin 
magnitude \a\ = 0.515. These configurations are summa- 
rized in Table [H 

Before discussing these configurations, we describe 
our numerical measurements which enable a quantita- 
tive comparison between the evolutions. First, we com- 
pare the luminosities as functions of gravitational wave 
frequency, as this is an observable and allows for a direct 
comparison of the different cases. We obtain frequencies 
from the I = 2, m = 2 gravitational mode. 

Second, we compute three different luminosities for 
each case: (i) the collimated luminosity ^collimated ob- 
tained by integrating the electromagnetic flux over a cone 
of points within 15° from the center of mass of the sys- 
tem; (ii) the non-collimated, or "isotropic," luminosity 
^isotropic obtained from the integral over an encompass- 
ing sphere minus the collimated luminosity of (i); (iii) the 
gravitational wave luminosity, Lqw- These different lu- 
minosities are displayed for the three binary configura- 
tions in Fig. [3] 

Consider first the 0/0 and u/d configurations which 
have essentially the same total angular momentum. Ex- 
tensive numerical simulations (see for instance [37]) and 
simple estimates [38 indicate that both binaries will 
merge into a final black hole with essentially the same 
spin (a ~ 0.67), and, so for late times after- merger, the 
expected jet structure should be quite similar, dictated 
by the standard BZ mechanism. 

The binary with individual spins aligned reaches a 
higher orbital velocity before merger than the previous 
two cases; thus, the expected maximum power should 
be higher. Moreover, the resulting final black hole spins 
faster (a ~ 0.8) and thus its BZ associated power will 
be higher than that of the previous cases. Fig. [3] illus- 
trates this expected behavior by presenting the Poynting 
flux energy vs. gravitational wave frequency for the three 
cases. 

As evident from the figure, at low frequencies where 
the orbital dynamics are the same in all cases, both spin- 
ning cases have a higher output than the non-spinning 
one and the difference is provided by the spin contri- 
bution to the jet emission. Furthermore, both spinning 
cases have equal collimated power output because the 
spin contribution to the luminosity depends only on the 
spin magnitude. 

We can further examine the basic relation explaining 
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the observed flux by comparing the three cases studied. 
We estimate the black hole coordinate velocities v and 
(collimated) luminosities at three representative frequen- 
cies Qi = {1;1.5;2}10 -5 [Hz/M 8 ] (i = 1..3) before the 
strong non-linear interaction starts. Since for these fre- 
quencies the measured speeds for the 0/0 and u/d cases 
are essentially the same, we notice that the flux of energy 
contributed by the spinning black holes can be estimated 
to be: 



L e Bz(ty — L u /d(P) - Z>o/o(^) 5 



(5) 



thus, 
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(A)/o; L u / d )(MsB4) 2 L 43 


L% st z (M 8 B 4 ) 2 L 43 
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(1.08; 1.62) 


0.54 


2 


(1.83; 2.30) 


0.47 


3 


(2.30; 2.80) 


0.50 



Notice that the remains fairly constant through 

these frequencies which is evident also in the Fig. [3] With 
this value one can estimate the luminosity for the u/u 

case as: 



+ Lf z (fl). (6) 



Applying such expression to our representative values 
we obtain, 
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(0.192; 0.174) 


(1.64; 1.85) 


2 


(0.206; 0.186) 


(2.52; 2.71) 
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(0.211; 0.187) 


(3.26; 3.42) 



thus the estimates are within ~ 10% of the measured 
values while the black holes are sufficiently separated that 
their jets do not strongly interact. Notice that at higher 
frequencies the aligned (u/u) case indeed has a higher 
associated power. 

In all cases, a significant non-collimated emission is in- 
duced (illustrated in the top right plot of Fig. [3| during 
the merger phase. Clearly, the simple minded picture of a 
jet produced by the superposition of the orbital and spin- 
ning effect can not fully capture the complete behavior 
at the merger epoch, although it serves to understand 
the main qualitative features and provides a means to 
estimate the power of the electromagnetic emission. 



IV. FINAL COMMENTS 

We have studied the impact of black hole motion 
through a plasma and indicated how the interaction can 
induce powerful electromagnetic emissions even for non- 
spinning black holes. Despite having examined a very 
small subset of the binary black hole parameter space, 
the results presented both here and in JT0j [11] suggest a 
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FIG. 3: Luminosity as a function of the orbital frequency 
for the binary black hole configurations described in Table [I] 
assuming M = 1O 8 M and B = 10 4 G. Top left: the colli- 
mated luminosity associated with the jets. Top right: the 
isotropic luminosity representing electromagnetic flux not as- 
sociated with the jets. Bottom: the gravitational wave out- 
put. 



wide applicability to general black hole binaries. More- 
over, as the plasma generally has a negligible effect on the 
dynamics of the black holes, then one needs only to know 
the dynamics of the black holes, say by numerical solu- 
tion or other approximate methods, in order to estimate 
the expected electromagnetic luminosity. A recent exam- 
ple is the work of [39] which uses the expected ~ B 2 v 2 
scaling for non-spinning black hole binaries coupled with 
the known distance dependence in time to obtain excel- 
lent agreement with the obtained luminosity from such a 
system. 

Naturally, spinning black holes produce stronger jets, 
and these jets, as shown earlier [TT] . will be aligned with 
the asymptotic magnetic field direct ion [42]. Thus, the re- 
sulting electromagnetic luminosity can be estimated to be 
Lfjj 4 ~ L S pin + £ S peed, which can be significant and have 
associated time- variabilities tied to the dynamical behav- 
ior of the system. In particular, the luminosity tied to 
the motion can result significantly higher than that tied 
to the spin. In addition, eccentric orbits and spin-orbit 
interactions driving orbital plane precessions can induce 
important variabilities that can aid in the detection of 
these systems. Furthermore, a significant pulse of nearly 
isotropic radiation is emitted during merger, thereby al- 
lowing observations of the system along directions not 
aligned with the jet. 

Consequently, binary black hole interactions with sur- 
rounding plasmas can yield powerful electromagnetic 
outputs and allow for observing these systems through 
both gravitational and electromagnetic radiation. Grav- 
itational waves from these systems corresponding to the 
last year before the merger could be observed to large 
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distances with LISA (up to redshifts of 5-10 [40 J) for 
masses ~ 1O 4_T M0 or earlier in the orbiting phase 
(and possibly through merger) via Pulsar Timing Ar- 
ray observations [41], targeting binaries with masses in 
M ~ 1O T_1O M0. As we have indicated here, both scenar- 
ios can have strong associated electromagnetic emissions. 
Our inferred luminosities of several 10 43 (M8l?4) 2 ergs/s 
corresponds to an isotropic bolometric flux of F x 10 -15 
erg (M 8j B 4 ) 2 /(cm 2 s) that could be detected to redshifts 
of z « 1 and even further depending on anisotropics, de- 
pending on the efficiency of processes taping this avail- 
able energy and producing observable signals. 



FIG. 4: Electromagnetic flux corresponding to the 0/0 (top 
row) and u/d (bottom row) cases at times ( — 16.7,-2.5) 
hrsMg and (—10.5,3.6) hrsMs (with respect to the merger 
time as marked by the peak in strength of gravitational wave 
emission) respectively. The orbiting stage leaves its imprint 
as twisted tubes. 
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FIG. 5: Electromagnetic flux corresponding to the u/u case at 
times (—2.5,11.6,40,68) firs Ms (with respect to the merger 
time as marked by the peak in strength of gravitational wave 
emission). The orbiting behavior leaves a clear imprint in 
the resulting flux of energy and, as the merger proceeds emis- 
sion along all directions is evident. At late times, the system 
settles to a single jet as dictated by the BZ mechanism. 
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